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Abstract 

Modular and quasimodular forms have played an important role in gravity and string the- 
ory. Eisenstein series have appeared systematically in the determination of spectrums and 
partition functions, in the description of non-perturbative effects, in higher-order corrections 
of scalar-field spaces, . . . The latter often appear as gravitational instantons i.e. as special 
solutions of Einstein's equations. In the present lecture notes we present a class of such 
solutions in four dimensions, obtained by requiring (conformal) self-duality and Bianchi 
IX homogeneity. In this case, a vast range of configurations exist, which exhibit interest- 
ing modular properties. Examples of other Einstein spaces, without Bianchi IX symmetry, 
but with similar features are also given. Finally we discuss the emergence and the role of 
Eisenstein series in the framework of field and string theory perturbative expansions, and 
motivate the need for unravelling novel modular structures. 



To appear in the proceedings of the Besse Summer School on Quasimodular Forms - 2010. 
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1 Introduction 

Modular forms often appear in physics as a consequence of duality properties. This comes 
either as an invariance of a theory or as a relationship among two different theories, under 
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some discrete transformation of the parameters. The latter transformation can be a simple 
Z2 involution or an element of some larger group like SL (2, Z) . The examples are numerous 
and have led to important developments in statistical mechanics, field theory, gravity or 
strings. 

One of the very first examples, encountered in the 19th century, is the electrtic-magnetic 
duality in vacuum Maxwell's equations (see e.g. [1]), which are invariant under interchang- 
ing electric and magnetic fields. This was revived in the more general framework of Abelian 
gauge theories by Montonen and Olive in 1977 |2j and culminated in the Seiberg-Witten du- 
ality in supersymmetric non- Abelian gauge theories |3j. There the duality group SL(2,Z) 
acts on a complex parameter r = S- + ^jr, where 8 is the vacuum angle and g is the coupling 
constant. The modular transformations give thus access to the non-perturbative regime of 
the field theory. 

In statistical mechanics, the Kramers-Wannier duality [4[ predicted in 1941 the existence 
of a critical temperature T c separating the ferromagnetic (T < T c ) and the paramagnetic 
(T > T c ) phases in the two-dimensional Ising model. The canonical partition function of 
the model, computed a few years later by Lars Onsager |5j, is indeed expressed in terms of 
modular forms. 

Over the last 30 years, modular and quasimodular forms have mostly emerged in the 
framework of gravity and string theory. At the first place, one finds (see e.g. |6j) the canonical 
partition function of a string of fundamental frequency co at temperature T: 

Z = (1.1) 

where q = exp 2inz = exp —h^/kT and tj(q) the Dedekind function (we refer to the appendix 
for definitions and conventions on theta functions, ( A.l |-( A.5| l). The average energy stored 
in a string at temperature T is thus 

9 In Z fico „ , , .„ „. 

< E > " " aiAF " II E2(Z) ' (1 ' 2) 

where £2(2) the weight-two quasimodular form. Again, the modular properties of these 
functions translate into a low-temperature/high-temperature duality, which exhibits a criti- 
cal temperature, signature of the Hagedorn transition. 

In the above examples the modular group acts on a modular parameter related to the 
temperature. There is a plethora of examples in gravity and string theory of more geometri- 
cal nature, related to gravitational configurations and in particular to instantons. 

An instanton is a solution of non-linear field equations resulting from an imaginary-time 
i.e. Euclidean action S [<p] as 
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It should not to be confused with a soliton. The latter is a finite-energy solution of non-linear 
real-time equations of motion and appear in a large palette of phenomena such as the prop- 
agation of solitary waves in liquid media (as e.g. tsunamis^} or black holes in gravitational 
set ups. 

Instantons have finite action and enter the description of quantum-mechanical processes, 
which are not captured by perturbative expansions, as their magnitude is controlled by 
exp — Vg at small coupling g (in electrodynamics g = e2 /hc). These phenomena include 
quantum-mechanical tunneling and, more generally, decay and creation of bound states. 
Their amplitude is weighted by exp(— s /h), where S is the action of the instanton solution 
interpolating between initial and final configurations (see (8| for a pedagogical presentation 
of these methods). 

All interacting (i.e. non-linear) field theories exhibit instantons. These emerged origi- 
nally in Yang-Mills theories (9 [10 [ as well as in general relativity |TT|-[13"|. In the latter case, 



their usefulness for the description of quantum transitions is tempered by quantum incon- 
sistencies of general relativity. Such configurations turn out nevertheless to be instrumental 
in modern theories of gravity, supergravity and strings for at least two reasons. 

At the first place, some gravitational instantons falling in the class of asymptotically lo- 
cally Euclidean (ALE) spaces have the required properties to serve as compactification set 
ups for superstring models usually defined in space-time of dimensions 10. This is the 
case, for example, of the Eguchi-Hanson gravitational instanton |12 13j, which appears as 
a blow-up of the C2/Z2 Ai-type singularity, or of more general Gibbons-Hawking multi- 
instantons [14| . 

The second reason is that supergravity and string theories contain many scalar fields 
called moduli. Their dynamics is often encapsulated in non-linear sigma models, which 
happen to have as a target space certain gravitational instantons such as the Taub-NUT, 
Atiyah-Hitchin, Fubini-Study Pedersen or Calderbank-Pedersen spaces JTT 15-21 1. Due to 
some remarkable underlying duality properties, most of the spaces at hand are expressed in 
terms of (quasi) modular forms, and this makes them relevant in the present context. 

It should be finally stressed that in the framework of string and supergravity theories, 
quasimodular forms do not appear exclusively via compactification or moduli spaces. Re- 
cent developments on the perturbative expansions in quantum field theory reveal how rele- 
vant the spaces of quasimodular forms are for understanding the ultraviolet behaviour and 
its connections with string theory acting as a ultraviolet regulator [22|. They also call for 
introducing new objects, which stand beyond the realm of Eisenstein series |23|. 



J A valuable account of these properties can be found in \ 7\. 
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2 Solving Einstein's equations 



It is a hard task in general to solve Einstein's equations. In four dimensions with Euclidean 
signature, following the paradigm of Yang-Mills, the requirement of self-duality (or of a 
conformal variation of it) often leads to integrable equations. Those are in most cases re- 



lated to self-dual Yang-Mills reductions, and possess remarkable solutions (see. e.g. |24|). It 
should be mentioned for completeness that self-duality can also serve as a tool in more than 
four dimensions. In seven or eight dimensions, it can be implemented using G2 or quater- 
nionic algebras |25-28j. It is not clear, at present, whether in those cases some interesting 
and non-trivial relationship with quasimodular forms emerge. We will therefore not pursue 
this direction here. 

2.1 Curvature decomposition in four dimensions 

The Cahen-Debever-Defrise decomposition, more commonly known as Atiyah-Hitchin- 
Singer 1 29 30j, is a convenient taming of the 20 independent components of the Riemann 



tensor. In Cartan's formalism, these are captured by a set of curvature two-forms (a,b, . . . = 
1,...,4) 

K\ = dw a b + < A cv\ = l -R\ cd d c A 6 d , (2.1) 

where {9 a } are a basis of the cotangent space and to a b = T a b 6 C the set of connection one- 
forms obeying the requirement of vanishing torsion 

r a = de a + w a b a e b = -T% c e b a e c = o. (2.2) 

The cyclic and Bianchi identities (d A d6 a = d A dco a b = 0), assuming a torsionless connec- 
tion, read: 

K\NB h = 0, (2.3) 
dVf b + co a c A 1Z c b — lZ a c A co c b = 0. (2.4) 

We will assume the basis {9 a } to be orthonormal with respect to the metric g 

g = S ah d a d h , (2.5) 

and the connection to be metric (Vg = 0), which is equivalent to 

co ah = -w ba , (2.6) 



The latter together with ( |2.2[ | determine the connection. 

The general holonomy group in four dimensions is SO (4), and g is invariant under local 
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transfromations A(x) such that 

m' — a -l' 

b 

under which the connection and curvature forms transform as 



a ' = A- la h 9 h , (2.7) 



< = A- la c w c d A\ + A- la c dA c b/ (2.8) 
n = A- la c K c d A d b . (2.9) 

Both o; fl /, and 1Z a \, are antisymmetric-matrix-valued one-forms, belonging to the representa- 
tion 6 of SO (4). 

Four dimensions is a special case as SO(4) is factorized into SO(3) x SO(3). Both con- 
nection and curvature forms are therefore reduced with respect to each SO (3) factor as 
3x1 + 1x3, where 3 and 1 are respectively the vector and singlet representations (z, ;',...= 
1,2,3): 

£i = \ (tooi + ^e ijk wA , A; = - f^oi - -e^a^ , (2-10) 
5 f = ^ fax + ^fc^) , ^i = \ fai - \e ijk nA , (2.11) 
while ( |2.1| | reads: 

5,- = dE; - e yfc E' A Y}, Ai = dAj + e ijk A' A A k . (2.12) 

Usually (Ej, 5,) and (A,-, A.,) are referred to as self-dual and anti-self-dual components of 
the connection and Riemann curvature. This follows from the definition of the dual forms 
(supported by the fully antisymmetric symbol e„f, C( Q 

&\ = \e\c d " C d. (2-13) 
= ^V^V (2-14) 

borrowed from the Yang-Milla^J Under this involutive operation, (E ; , <S ; ) remain unaltered 
whereas (A,, Ai) change sign. 



2 A remark is in order here for D = 7 and 8. The octonionic structure constants t/> a fi 7 a, f>, J G {1/ • • • , 7} and 
the dual G2-invariant antisymmetric symbol xp K ^ s allow to define a duality relation in 7 and 8 dimensions with 
respect to an SO(7) D G2, and an SO(8) D Spin 7 respectively. Note, however, that neither SO(7) nor SO(8) is 
factorized, as opposed to SO (4). 

3 Note the action of the duality on the components, as cv" b = T\ c 8 c , H\ = \R a bcd c A 6 d : 

pa 1„* fre 

1 he ~ 2 fc he 1 fe> 



be 2 be fc 

b 

and similarly for the Weyl part or the Riemann. 



R" — If" ^R e 

R bed ~ 2 fc be lx fed' 
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Following the previous reduction pattern, the basis of 6 independent two-forms can be 
decomposed in terms of two sets of singlets /vectors with respect to the two SO (3) factors: 

<p = 6° A &■ + hjpfii A 6 k , (2.15) 

x l = e° AP-~j }k ei A6 k . (2.16) 

In this basis, the 6 curvature two-forms 5, and A[ are decomposed as 



S\ = r 1? 

A 2 U 



(2.17) 



where the 6x6 matrix r reads: 



A C+\ (W+ C+\ s , v 

r=\ = + 7 1 6 - 2.18 

\c- By yc- w-y 6 

The 20 independent components of the Riemann tensor are stored inside the symmetric 
matrix r as follows: 

• s = Tr r = 2Tr A = 2Tr B = R /2 is the scalar curvature. 

• The 9 components of the traceless part of the Ricci tensor S a b = R a b — jg a b (Rab = ^ c flC i,) 
are given in C + = (C - )' as 

S 00 = TrC+, S oi = e i ik C7 k , S tj = C± + C£ -TrC+fy (2.19) 

• The 5 entries of the symmetric and traceless W + are the components of the self-dual 
Weyl tensor, while W~ provides the corresponding 5 anti-self-dual ones. 

In summary, 

S l = Wf + tflH + lc+x', (2-20) 

Ai = Wr + ^ + ^C^, (2.21) 

where 

W? = \W$¥. Wr = ^W^ (2-22) 

are the self-dual and anti-self-dual Weyl two-forms respectively. 

Given the above decomposition, some remarkable geometries emerge (see e.g. [31] for 
details): 

Einstein C± = (o R ah = \g ai ) 
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Ricci flat C* = 0, s = 

Self-dual A, = O {W~ = 0, C± = 0, s = 0} 

Anti-self-dual 5 i = 0«{W+=0,C ± =0,s = 0} 

Conformally self-dual W = 

Conformally anti-self-dual W + = 

Conformally flat W+ = W = 

Note that self-dual and anti-self-dual geometries are called half-flat in the mathematical lit- 
erature, whereas self-dual and anti-self-dual is meant to be conformally self-dual and anti- 
self-dual. 



2.2 Einstein spaces 

The self-dual and anti-self-dual geometries have a special status as they are automatically 
Ricci flat: 

Ai = or Si = => C ± = 0, s = 0. (2.23) 

They provide therefore special solutions of vacuum Einstein's equations, which include 
gravitational instantons already quoted in the introduction such as Eguchi-Hanson, Taub- 
NUT or Atiyah-Hitchin. 

More general solutions are obtained by demanding conformal self-duality on Einstein 
spaces 

W+ = or W~ = and C ± = (2.24) 
with non-vanishing scalar curvature^] 

s = 2A. (2.25) 

Those are the quaternionic spaces and include other remarkable instantons such as Fubini- 
Study Pedersen or Calderbank-Pedersen. 

Conditions ( 2.24} and ( 2.25} can be elegantly implemented by introducing the on-shell 
Weyl tensor 

W ah = Tl ab - j6 a A e h . (2.26) 
These 6 two-forms can be decomposed into self-dual and anti-self-dual parts: 

W+ = 5 i -^ = W+ + ^(s-2A)^ + ic+^ / (2.27) 
% = A-^ = Wr + l(s-2Ak, + ^C ! 7^. (2.28) 



4 Requiring vanishing amounts to demanding the space to be Einstein (R a f, = ^g„b), which implies that 
its scalar curvature is constant. 
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Quaternionic spaces are therefore obtained by demanding 



W+ = or Wr = 0. 



(2.29) 



Furthermore, using the on-shell Weyl tensor ( 2.261 , the Einstein-Hilbert action reads: 



Seh = 



32nG Jm 



e a bcd [ W ab + A 6° ) A 6 C A 



A 



(2.30) 



3 Self -dual gravitational instantons in Bianchi IX 



Inspired by applications to homogeneous cosmology (see e.g. |32 [), spaces M.± topologically 
equivalent to R x have been investigated extensively in the cases where M3 are homo- 
geneous of Bianchi type. These foliations admit a three-dimensional group of motions acting 
transitively on the leaves 

The study of all Bianchi classes (I-IX) has been performed (for vanishing cosmological 
constant) in 1 33 -[35 [ and more completed recently in p6| . It turns out that only Bianchi IX 
exhibits a relationship with quasimodular forms. 

3.1 Bianchi IX foliations 

Under the above assumptions, a metric on M.^ can always be chosen as (see e.g. (37j) 

ds 2 = dt 2 + gij {tyai, (3.1) 

where a 1 , i = 1, 2, 3 are the left-invariant Maurer-Cartan forms of the Bianchi group, satisfy- 
ing 

do* = y jk a'Aa k . (3.2) 
This geometry admits three independent Killing vectors tangent to Ai$ and such that 

[ZuSj] =c i ] £ k . (3.3) 

In the case of Bianchi IX, the group is SU(2). Using Euler angles, the Maurer-Cartan 
forms read: 

a 1 = sin d sin ip dip + cos ip dd 

a 2 = sin # cos t/; dtp — sini/?d# (3.4) 
o 3 = cost? dip + dtp 



with 0<#<7r,0<<p< 2n,0 < ip < 47T. The structure constants are c ! - 



(jk 
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with £123 = 1. Similarly the Killing vectors are 



£ x = - sin (p cot #3^ + cos(p3# + ^|9^ 
£2 = cos <p cot d dtp + sin <p d$ - 3^ 
£3 = dtp- 



OS) 



Although for some Bianchi groups it is necessary to keep in ( |3.1| general, for Bianchi 
IX it is always possible to bring it into a diagonal form, without loosing generality (for a 



systematic analysis of this, see |36]). We will make this assumption here, introduce three 
arbitrary functions of time Q ! as well as a new time coordinate defined as dt = V Q 1 Q 2 Q 3 dT, 
and write the most general metric ( |3.1| on a Bianchi IX foliation as 



Q 1 V J Q2 
For this metric, the two-form basis ( 2.15 1 and ( 2.16| l reads: 



Q3 



CVCi k dT A (f + A a k , 

an k dT a cr» - qv a ^, 



(3.6) 



(3.7) 
(3.8) 



where z, /, A; are a cyclic permutation of 1, 2, 3 without over i. Using Eqs. ( |2.2| , ( |2.6| and ( 2.10| , 
one finds for the corresponding Levi-Civita connection 



Ei = 
A, = 



4v / n 1 Q 2 n 3 V n* - ~ n/ ~ o fc 
1 /fi''-fv'fi fc (V'-nW tf-nw" 



4v / n 1 Q 2 Q 3 



(3.9) 
(3.10) 



where / stands for d //dr (as previously, i,j,k are a cyclic permutation of 1,2,3 and no sum 
over i is assumed) . 

3.2 First-order self-duality equations 

From now on, will focus on self-dual solutions of Einstein vacuum equations (anti-self-dual 
solutions are related to the latter e.g. by time reversal). Following ( 2.10| , self-duality equa- 



tions (2.231 read: 



dAi + e ijk A! A A k = 0. 



(3.11) 



Equations (3.11 1 are second-order. They admit a first integral, algebraic in the anti-self-dual 



connection A;: 



A{ = -^-cri with Aij = or 5{y 



(3.12) 
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Put differently, vanishing anti-self-dual Levi-Civita curvature can be realized either with a 
vanishing anti-self-dual connection, or with a specific non-vanishing one that can be set to 
zero upon appropriate local SO (3) C SO(4) frame transformation (see J3l) for a general 
discussion, [38 1 for Bianchi IX, or |36| for a more recent general Bianchi analysis). These two 
possibilities lead to two distinct sets of first-order equations. In the present case, using ( 3.10} 
one obtains: 

A,- = 0o jo 1 = Q 2 Q 3 , Q 2 = Q 3 Q 1 , Q 3 = Q 1 Q 2 | , (3.13) 



and 



Aj = Sjj— < 



2 



Q 1 = Q 2 Q 3 -Q 1 (Q 2 + Q 3 ) 

Q 2 = Q 3 Q 1 - Q 2 (Q 3 + Q 1 ) (3.14) 
Q 3 = Q 1 Q 2 -Q 3 fQ 1 + Q 2 ) . 



Historically, both systems were studied in the 19th century in the search of integrals lines 
of vector fields. The first is the Lagrange system, appearing as an extension of the rigid- 
body equations of motion. It is algebraically integrable and was solved a la Jacobi. The 
second set is called Darboux-Halphen and appeared in Darboux's work on triply orthogonal 
surfaces |39| . Generically, it does not possess any polynomial first integral, and was solved 
by Halphen in full generality using Jacobi theta functions [40 1. 

In the late seventies, integrable systems of equations such as Lagrange or Darboux- 
Halphen emerged in a systematic manner in self-dual Yang-Mills reductions |24| . This has 
led many authors to investigate these equations in great detail and, in particular, to unravel 
their rich integrability properties (see e.g. |41 -43j as a sample of the dedicated literature). It 



took a long time, however, to realize that these systems were actually related with gravita- 
tional instantons, foliated by squashed spheres. 

When all three Q's are identical, the leaves of the foliation are isotropic three-spheres 
with SU(2) x SU(2) isometry generated by the above left Killing vectors = 1,2,3 ( |3.5| , 
as well as by three right Killing vectors 

( e\ = - sin ip cot d + cos xp d# + 3<p 
e 2 = -cos xp cot tidy -sinxpd# + d 9 (3.15) 
63 = dtp- 

Lagrange and Darboux-Halphen systems are equivalent in this case (actually related by time 
reversal), Q = ±Q 2 , and the solutions lead to flat Euclidean four-dimensional space. 

More interesting is the case where Q 1 = Q 2 ^ Q 3 . Here, the leaves are axisymmet- 
ric squashed three-spheres, invariant under an SU(2) x IT(1) isometry group generated by 
i = 1, 2, 3 and 63. On the one hand, the Lagrange system leads to two distinct gravitational 
instantons known as Eguchi-Hanson I and II ]12|13| , out of which the first has a naked singu- 
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larity and is usually discarded. On the other hand, the Darboux-Halphen equations deliver 
the celebrated Taub-NUT instanton (TTJ. 

Thanks to the algebraic integrability properties of Lagrange system, it took only a few 
months to Belinski et al. to generalize the Eguchi-Hanson solution to the case where Q 1 ^ 



Q 2 Q 3 [44] - the symmetry is strictly SIT (2) but the solution is plagued with naked singu- 



larities. A similar generalization of the Taub-NUT solution turned out much more intricate, 



and after some fruitless attempts |38[, Atiyah and Hitchin reached a regular solution, eli- 
gible as a gravitational instanton and expressed in terms of elliptic functions [15[. It was 
only realized in 1992 by Takhtajan [41 [ that first-order self-duality equations for Bianchi IX 
gravitational instantons were in fact Lagrange and Darboux-Halphen systems, and that the 
Atiyah-Hitchin instanton was a particular case of the general solution found by Halphen in 
1881 (40). 

It is finally worth mentioning that the above systems of ordinary differential equations 
also appear in the framework of geometric flows in three-dimensional Bianchi IX homoge- 
neous spaces. The original mention on that matter can be found in [45 ]; later and indepen- 
dently it was also quoted in |46) . At that original stage, this relationship was limited to the 
case of Bianchi IX with diagonal metric. It was proven recently to hold in full generality in 
all Bianchi classes ]47][48} . 



4 The Darboux-Halphen system 

The Darboux-Halphen branch of the self-duality first-order equations of Bianchi IX folia- 
tions in vacuum is the most interesting for our present purpose as it is the one related with 
quasimodular forms. 

4.1 Solutions and action of SL(2, C) 

Consider the system in the complex plane: oo l {z), z G C satisfying 



dw 1 
dz 

dw 2 
dz 

dw 3 
dz 



oo 2 oo 3 — GO 1 (go 2 + GO 3 ) 

oo 3 co l - cv 2 (go 3 + a; 1 ) (4.1) 

G0 1 G0 2 — CO 3 (cO 1 + GO 2 ) . 



The general solutions of this system have the following properties |40 41}: 

• The <x>s are regular, univalued and holomorphic in a region with movable boundary 
{i.e. a dense set of essential singularities). The location of this boundary accurately 
determines the solution. 
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If co l {z) is a solution, thus 



Q\z) = ~^— 2 cv> ( + a MgSL(2 /C ) (4.2) 

(cz + df \cz + dj cz + d' \ c dj V ' 



is another solution^] with singularity boundary moved according to z 



az+b 
cz+d ' 



The resolution of the equations and the nature of the solutions strongly depend on whether 
the cos are different or not. In the case where co 1 = co 2 = oo 3 , the solution is simply 

„W = _J_ (4.3) 

Z-Zq 



with Zo an arbitrary constant. Under SL(2,C), the new solution o)'(z) is of the form (4.3 
with the pole displaced according to 

dzo -b ,. 
z = ■ 4.4 

czq — a 

If oo l = co 2 co 3 the solutions are still algebraic: 

/ ,1/2 1 , ,3 Z ~ Z * (A r\ 

Z-Zq [Z-ZqY 

with two arbitrary constants: zo,z*. A simple pole for co 1,2 , and double for oo 3 appears at 
zo, whereas z* is a root for to 3 . Acting with SL(2,C) keeps the structure \A.5\ with new 
parameters: 

dzo-b _ _ zo-z* 

z = , z - z t . = -j. (4.6) 

czq - a ( CZq _ a f 

The fully anisotropic case is our main motivation here. In this case no algebraic first 
integrals exist and the general solution (see |40-42j) is expressed in terms of quasimodular 
forms, oo l S QM\ (r(2)), where T(2) is the level-2 congruence subgroup of SL(2,Z) (the 
subset of elements of the form ( a c b d ) = (o?) mod 2) . Concretely 

a/(z) = -^log^( Z ) (4.7) 



5 The same property holds for the Lagrange system {3.13 \, limited to the subgroup of transformations of the 

a b 

o v. , 



form ( |] i/„ )• This solution-generating pattern based on the SL(2,R) is closely related to the Geroch method 



|49]|50 
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with £ t (z) triple^ of holomorphic zveight-2 modular forms of T(2). Again, the SL(2,C) ac- 



tion (4.2 1 generates new solutions {co*} — > {&'} with a displaced set of singularities in C, 
whereas the SL(2,Z) C SL(2,C) acts as a permutation on ws. 

Note for completeness that real solutions of the real coordinate T are obtained from the 
general solutions as 

Q^(T) = icv e (iT) = -1 Alegar)- (4-8) 



According to (4.2 1, new real solutions are generated as 



4.2 Relationship with Schwartz's and Chazy's equations 

Anisotropic solutions of the Darboux-Halphen system (a; 1 ^ co 2 ^ co 3 ) exhibit relationships 
with other remarkable equations. Define 

A = (4.10) 

CO 1 — CO 

For to 1 solving Darboux-Halphen equations, A is a solution of of Schwartz's equation 
A'" 3 fk"\ 2 1/1 1 1 \ /w , 2 



Conversely, any solution of the latter equation provides a solution for the Darboux-Halphen 
system as the following triplet: 

-l_«»Vdz 2 _dA/dz c3 _ d V dz / d1 ,x 

* " A ' fc ~A-1' * ~A(A-1)' (4 ' iZj 
from which it is straightforward to show that 

£ 1 - £ 2 + £ 3 = 0. (4.13) 

We also quote for completeness 



co 1 — co 2 1 — A co 3 — to 2 



1 — A co 3 — co 2 ' A a; 1 — a; 3 ' 
6 Notice their general transformations as generated by z — > — l/z and z + 1 : 

2 -> -Vz : (fii £jj £») -> z 2 £"l -£3) 

z ->■ z + 1 : £" 2 £3) - (£3 £" 2 £1) • 



(4.14) 
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Define now 

y = -2 (to 1 + to 2 + to 3 ^ . (4.15) 
Again, for solutions of Darboux-Halphen equations to', y a solution of Chazy's equation p5l) 

j/"=2yy"-3(y') 2 . (4-16) 



The first and second derivatives of y provide the remaining symmetric products 



y' = 2 (co 1 ^ 2 + to 2 to 3 + to 3 co^ , 



" - -\2co 1 co 2 to 3 . 



The Jacobian relating {co 1 , co 2 , co 3 } to {y,y',y"}, 

J =(to 1 - co 2 )(to 2 - to 3 )(to 3 - co 1 ) 



(4.17) 
(4.18) 



(4.19) 



is regular for co 1 ^ co 2 ^ co 3 . The latter are alternatively obtained by solving the cubic 
equation 



oo 3 + ^yto 2 + ^y'co + ^y" = 0, 



(4.20) 



for any solution y of Chazy's equation. 



4.3 The original Halphen solution 



A particular solution of the Darboux system ( 3.14| is the original Halphen solution |40| . In 
this language, it corresponds to Ah = ^/ 



S 1 
E 2 
E 3 



ind\ 
—in-d\ 



I (E 2 + ^ - 



(4.21) 



CO' 



H 



This is also the solution found by Atiyah and Hitchin [151 as the Bianchi IX gravitational 
instanton solution relevant for describing the configuration space of two slowly moving 



BPS SU(2) Yang-Mills-Higgs monopoles |52 53J. The corresponding Chazy's solution and 
derivatives are combinations of (holomorphic) Eisenstein series (see appendix, Eqs. ( A.6| ): 



ircEx 

E 4 ) 

% (E 3 -3E 2 E^ + 2E 6 ). 



in (p2 
6 1^2 

T 3 



(4.22) 



Starting from ( 4.21| | all solutions are obtained by SL(2,C) action (4.2 
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5 Back to Bianchi IX self-dual solutions 

Any real solution {Q'(T)} of the Darboux-Halphen system provides a four-dimensional 
self-dual solution of Einstein vacuum equations in the form ( |3.6| . Not all these solutions are 
however bona fide gravitational instantons as some regularity requirements must be fulfilled. 



5.1 Some general properties 

An elementary consistency requirement is that the metric ( |3.6| should not change sign along 
T. In particular, a simple root of a single Q ! turns out to be a genuine curvature singularity. 
Assuming e.g. linearly vanishing Q 1 and introducing as time coordinate the proper time T 
around the root, the metric locally reads: 

ds 2 « dr 2 + A (V) 2 + Yr 2/3 ( (a 2 ) 2 + (a 3 ) 2 ) (5.1) 



T 



with S, Y constants. This metric has a curvature singularity at T = 0. 

Other pathologies can appear, which do not necessarily affect the consistency of the solu- 
tion. Poles of some Qs or multiple roots are potential natural boundaries or (non-)removable 
coordinate singularities such as bolts or nuts. The latter are fixed points of some Killing vec- 
tors £ (V( v £^) = 0), for which the matrix V^c^j is respectively of rank 2 and 4. A general, 
complete and detailed presentation of these properties is beyond the scope of these notes 
and is available in the original paper |54| . For our purpose here, we recall two generic sit- 
uations, where again we present the metric in local proper time T around a fixed point at 
r = 0: 



Rank 2 - bolt This singularity is removable if the metric behaves as 



ds 2 « dT 2 + C 2 ((^) 2 +(<7 2 ) 2 )+^(^ 



x 2 



= dr 2 + ^f- (dip + cos ddcp) 2 + I 2 (d# 2 + sin 2 ddcp 2 ) 

and provided nx P/i S [0, 2n[. Locally the geometry is thus R 2 x S 2 . 
Rank 4 - nut This singularity is removable if the metric behaves as 

ds 2 « dT2 + £ ((rf+irf+io*) 2 ) 

= dr 2 + £ ( di? 2 + sin 2 ddq> 2 + {dtp + cos dd(pf 
For a nut, the local geometry is R 4 (here in polar coordinates). 



(5.2) 



(5.3) 
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5.2 Behaviour of Darboux-Halphen solutions 

As already pointed out, there are thee distinct cases to consider: Q ! all equal, Q 1 = Q 2 ^ Q 3 
or Q 1 ^ Q 2 ^ Q 3 . In the first case, 

! ' = — ^— Vi, (5.4) 
1 — lo 

and the four-dimensional solution corresponds to flat space. When only two Qs are equal, 
the isometry group is extended to SU(2) x U(l) and real solutions read: 

Ql,2 = _L_ Q 3 = Z~l* . (5.5) 

There are 3 special points: T = T*,Tq and T — > oo. One can analyze their nature by zooming 
around them, using proper time: 

• At T — > oo one recovers the behaviour ( |5.3| and this point is a nut. 

• Around T = To one finds 

ds 2 ps dr 2 + t 2 (d<? 2 + sin 2 ddcp 1 ) + — (dip + cos Mcpf . (5.6) 

To — T* 

This is an S 1 fibration over R 3 , the fiber being dtp + cos ddcp. It is called Taubian infinity 
(see 1 54 [) and appears as a natural "boundary". 

• At T = T* there is a curvature singularity as the metric behaves like ( |5.1| with E = 

( ? \ 2/3 y — ( - y /3 

\3(T -T t )2) ' \2(T -T t )J ■ 

One therefore concludes that in order to avoid the presence of naked singularities, the singu- 
lar point T = T* should be hidden behind the Taubian infinity i.e. T* < To (see Fig. [TJ. Under 
this assumption, the self-dual solution at hand is well behaved and provides the Taub-NUT 
gravitational instanton 1 11 31). It is most commonly written as: 

~X + I ? - ft" 1 )' + ^) + (5-7) 



where m 2 = T \ T > and m(r — m) — 

The case where Q 1 ^ Q 2 7^ O 3 is the most interesting in the present context since it 
involves quasimodular forms. The real Halphen solution (see Eq. ( |4.21 1 with z = iT or 
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Figure 1: Generic solution Q 1 = Q 2 < O 3 . 



q = exp — 2nT)) reads: 



n^(T) = f (£ 2 + 1? 3 4 + 1? 4 4 ) 

o 3 I (r) = f (E 2 + #f-#f) <Q 2 1 . 



It is defined for T > with a pole at T = 0: 



L2 H~ 2T 2 



''< "2,3 _ 1 



Around this pole, the behaviour of the metric is 



ds 2 



dT 2 + T 2 ((^) 2 + ( £ , 2 ) 2 ) + ^(^) 2 ) 



(5.8) 



(5.9) 



(5.10) 



and we recover a Taubian infinity (S 1 fiber over R 3 ). The large-T behaviour is exponential 
towards a constant 



Q^ 3 « zp4 n eX p _ ttT, Qh ~ n /i + 4zr exp -2nT 



(5.11) 



with 



ds 2 



dr 2 + f ((^) 2 +(^) 2 )+4r 2 (^ 2 ) 2 ) 



(5.12) 
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This is precisely a bolt as in Eq. (5.2 1 with n = 4, £ = 
directions 2 and 3. All this is depicted in Fig. [2] 



\J n /i and permutation of principal 




Figure 2: Halphen original solution (Q^ < < < Q^). 

As already quoted, the self-dual vacuum geometry corresponding to Halphen's original 
solution is the Atiyah-Hitchin gravitational instanton [15|. Using modular transformations 



(4.9 \, one constructs all other real solutions with strict SU(2) isometry (i.e. with all Q ! differ- 



ent): 

n ' (T) ^?^ n «(^) + Cf^D' (5 ' 13) 

Are those well behaved? 

The answer is no because a root of one Q always appears between the Taubian infinity 
and the bolt. This root is a curvature singularity which spoils the regularity of the solution. 



In order to elaborate on that, we first observe that in Eq. (5.13 ', ct+d mus t be positive, as 



real are only defined for positive argument. Assume for concreteness that 



t^oo CT + D C 



On the one hand, at large T 



Q ! ' = ^+0(Vr2), (5.15) 

and trading T for the local proper time one finds a nut (see ( |5.3| ). On the other hand, the 
values Too = — D /c < Tq = — B / a correspond to two poles, and Q ! (T) are defined for T < T^ 
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or To < T with reflected behaviour. For To J$ T 

1 



n 1 



71 



2A 2 (T-T ) 2 ' 



n 



2,3 



r-r 



(5.16) 



(note the sign flip in Q 1 ) and 



ds 2 « dr 2 



(5.17) 



Therefore T = To is a Taubian infinity (S 1 fiber over R 3 ), and as T moves from T = To to 
T — > +oo one moves from the Taubian infinity ("boundary") to a nut. A similar conclusion 
is reached when scanning T from Too to — oo. 

The problem arises because there is always a value T* such that To < T* < oo with Q; = 
< Q 3 < Q 2 (see Fig. [3). This unavoidable root is a genuine curvature singularity of the 
metric. Because of this, no anisotropic solution of the Darboux-Halphen system other than 
the original one (( 4.21| or (5.8 1) provides a well-behaved Bianchi IX gravitational instanton. 




Figure 3: Generic solution for To < T and < O 1 < Q 3 < O 2 . 
5.3 A parenthesis on Ricci flows 

Ricci flows describe the evolution of a metric on a manifold, governed by the following first- 
order equation: 

= -% (5.18) 
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where Rh stands for the Ricci tensor of the Levi-Civita connection associated with gjj (see 
e.g. |55J). It was introduced by Hamilton in 1981 [56 1 in order to gain insight into the ge- 
ometrization conjecture of Thurston (see e.g. j57l), a generalization of Poincare's 1904 con- 
jecture for three-manifolds, finally demonstrated by Perel'man in 2003 |58|. Ricci flows are 
also important in modern physics as they describe the renormalization group evolution in 
two-dimensional sigma-models [59|. 

The case of homogeneous three-manifolds is important as it appears in the final stage 
of Thurston's geometrization. Homogeneous three-manifolds include all 9 Bianchi groups 
plus 3 coset spaces, which are H3, H2 x S 1 , S 2 x S 1 (S n and H n are spheres and hyperbolic 
spaces respectively) |60}6lJ. The general asymptotic behaviour was studied in detail in |62) . 



A remarkable and already quoted result |45 -48]| is the relationship between the parametric 
evolution of a metric 



2 /n 2 Q 3 / a 2 /Q3Q1 2 /QiQ 2 , 3,2 



ds = V IT V ) + V Hr M + V TP- (M9) 

on .M 3 of Bianchi type^J and the time evolution inside a self-dual gravitational instanton 
on M4 = R x A^3 as given in ( |3.6| : the equations are the same (f in ( 5.18} and T in ( |3.6| 



are related as dt = y Q 1 Q 2 Q 3 dT). Ricci flow on three-spheres is therefore governed by the 



Darboux-Halphen equations (3.14|. 



Solutions of the Darboux-Halphen system describe Ricci-flow evolution if Vz Q ! (T) > 0, 
assuming that this holds at some initial time To. It is straightforward to see that this is always 
guaranteed. Indeed, it is true when at least two Qs are equal, as one can see directly from 



the algebraic solutions (5.4! and (5.5 1. More generally, suppose that < Qj < Qq < Qq (the 



subscript refers to the initial time Tq) and that Qi has reached at time T\ the value Cl\ = 0, 



while Q 2 ,Q^ > 0. From Eqs. (3.14 1 we conclude that at time T%, C\\ = = —Cl\ Q3 < 
and Cl\ = Q2 > 0. This latter inequality implies that vanishes at T\ while it is 
increasing, passing therefore from negative to positive values. This could only happen if Qq 
were negative, which contradicts the original assumption. However, if indeed Qq < and 
Qq, Qq > 0, there is a time T\ where Q 1 becomes positive and remains positive together with 
Q 2 and Q 3 until they reach the asymptotic region. 

Solutions ( |5.4| and ( |5.5| show that the asymptotic behaviour of Qs is clearly V r, when at 
least two Qs are equal. In the more general case, the large-T behaviour is readily obtained 
thanks to the quasimodular properties of the solutions (see footnote [6j: 

n 1 ' 2 ' 3 (r) = -^Q 2 ' 1 ' 3 (i)+i. (5.20) 

7 The precise statement is actually formulated for more general, non-diagonal metrics, as explained in detail 
in 1 48 1, and is valid in all Bianchi classes. For Bianchi IX, the diagonal ansatz exhausts, however, all possibilities. 
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Therefore, for finite and positive Q = Q'(0), 



1 

Q' = — + subleading at large T, (5.21) 



as one observes in Fig. |4j Note that this does not hold for the solution (5.8 1 because for the 
latter T = is a pole and Qq = Q'(0) is neither finite, nor positive for all i. The behaviour at 
large T is not V T, but exponential (see Eq. ( 5.11} and Fig. [2j. 




Figure 4: Generic behaviour for < Qq < Qq < Qq. 



As a consequence of the generic behaviour ( 5.21| | of Qs for positive and finite initial con 



ditions, the late-time geometry on the S 3 under the Ricci flow is 

ds 2 «-^=((^) 2 +(0 2 +(- 3 ) 2 )- (5-22) 



This is an isotropic (round) three-sphere of shrinking radius^] It is worth stressing that this 
universal behaviour is specifically due to the quasimodular properties of the solution, re- 
flected in the non-co variant V T term of ( 5.20| . 



8 At large times, the original SU(2) or SU(2) X tl(l) isometry group gets enhanced to SU(2) x SU(2), while 
the volume shrinks to zero. These are generic properties along the Ricci flow: the isometry groups may grow in 
limiting situations, whereas the volume is never preserved, but shrinks for positive-curvature geometries: 
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6 Bianchi IX foliations and conformal self-duality 



So far, we have considered self-dual solutions of Einstein's equations. These satisfy Eqs. 



(2.231 and are Ricci flat. Solutions of the Darboux-Halphen system involving quasimodular 
forms are relevant in particular when Bianchi IX foliations are considered. The Lagrange and 
Darboux-Halphen systems, and more general modular and quasimodular forms emerge, 
however, in set ups where no self-duality and/or Bianchi IX foliation is assumed. Einstein 
conformally (anti-)self-dual spaces i.e. quaternionic spaces turn out to exhibit such interest- 
ing relationships. 

Conformally self-dual Einstein spaces satisfy (see Eqs. ( 2.29) ) 



wr = 0. 



(6.1) 



This two-form is defined in ( 2.28) as the anti-self-dual part of the on-shell Weyl tensor ( 2.26) . 



The latter includes a cosmological constant A and (6.1l implies that this space is Einstein 
(R a i, = Ag a t,) on top of being conformally self-dual (W~ = 0). 



6.1 Conformally self-dual Bianchi IX foliations 

Assuming the four-dimensional space be a foliation Ai^ = RxM3 with Ai^a general ho- 
mogeneous three-sphere invariant under SU(2) isometry, we can in general endow it with a 
metric ( |3.6) . The Levi-Civita connection one-forms of the latter are given in (3.9 1 and ( 3.10[ l). 
Conformal self-duality condition ( |6.1[ l does not require the flatness of the anti-self-dual com- 
ponent of the connection A, as in ( 3.11) . Hence, no first integral like ( 3.12) is available. 



In order to take advantage of the conformal self-duality condition (6.1 1 and reach first- 
order differential equations as in the case of pure self-duality, we can parameterize the con- 
nection A\ and E, and demand that (|6.1) be satisfied. This is usually done by setting both A, 



and E, proportional to a 1 , as in (3.12 \, with a T-dependent coefficient though (see e.g. [63 1): 



E, 



Bi i 
'ICY*' 

2Q' 



Using Eqs. (3.9 1 and (3.10 ', one obtains a relationship between {O'} and {B,, A,}: 



Q ! = Q'Q ,C - Q ! (A; + Ajt) = -tVn* + CT (By + B k ) . 



(6.2) 
(6.3) 

(6.4) 



Furthermore, Eqs. ( 2.27) and ( 2.28) lead to the following expressions for the on-shell Weyl 
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tensor: 



~ { 4QT^Q3 (* " BjB k - B t ( B} + B k )) + - A_ J tf, (6.5) 



+ { 4QWQ3 & " A, A, + A, (A, + A fc ) ) - £ } ^. (6.6) 

The additional (with respect to ( |6.4| |) first-order equations for {B,} or {A;} are obtained by 
imposing on-shell conformal self-duality. The canonical method for that is to demand that 



both coefficients of <p l and x' in (6.6 1 vanish. This guarantees (see ( 2.27 )) a conformally self- 
dual, Einstein manifold with scalar curvature R = 4 A, in other words a quaternionic space. 

Solving the system of equations obtained for conformally self-dual, Einstein manifolds 
depends drastically on whether or not the isometry is strictly SU(2), i.e. the leaves of the 
Bianchi IX foliation are anisotropic, triaxial spheres. When the isometry is extended to 
SIT (2) x U(l) (two equal Qs), the equations are algebraically integrable (as in the Darboux- 
Halphen system ( 3.14| ) and no relationship appears with modular or quasimodular forms. 



This leads to a variety of well known biaxial solutions (see p6|[17}[63| as well as |64| for a 
detailed presentation of the resolution) such as (anti-)de Sitter-Taub-NUT, (anti-)de Sitter- 
Eguchi-Hanson, (pseudo-)Fubini-Study - CP2, Pedersen (the parentheses correspond to 
negative A) . . . When all Qs are equal, the leaves are round, uniaxial three-spheres, and the 
only four-geometries are the symmetric S 4 or H4 (depending again on the sign of A). 

Although straightforward, the above approach leads for the triaxial case to equations 
which are not known to be integrable. Hence, their resolution is not systematic and general. 



An alternative strategy has been proposed by Tod and Hitchin [65.66|, based on twistor 



spaces and isomonodromic deformations (see also [67, 68 1). In a first step, one sets A to 



zero in (6.6 1 and demands the coefficient of x 1 to vanish. This is equivalent to demanding 
conformal self-duality and zero scalar curvature (W~ = s = 0) without setting C» to zero. 
Thus, the space is not Einstein and has zero scalar curvature. The final step is to perform a 
conformal transformation, which allows to restore a non-vanishing scalar curvature, while 
simultaneously setting =0. One thus obtains a quaternionic space. 

Explaining the details of this procedure is beyond our present scope, and we will there- 
fore limit our presentation to the issues involving modular forms, which stem out of condi- 
tions W~ = s = 0. These are imposed by demanding that the coefficient of x' vanishes in 
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( |6.6| > and setting A = 0: 
— / 

Ai = A 2 A 3 - Ai (A 2 + A3) 
I < A 2 = A 3 A 1 -A 2 (A 3 + A 1 ) 
^A 3 = AiA 2 -A 3 (Ax + A 2 ). 

They are supplemented with Eqs. (|6.4|> r which read for {A,}: 



(6.7) 



II 



Q 1 = Q 2 Q 3 -Q 1 (A 2 + A 3 ) 
Q 2 = Q 3 Q 1 - Q 2 (A 3 + Ai) 
Q 3 = Q 1 Q 2 -Q 3 (A1+A2). 



(6.8) 



Before pursuing the present investigation any further, it is worth making contact with the 



results of Sec. 3.2 on genuine self-duality equations. Assuming the system I and II satisfied 
i.e. W- = s = 0, Ai (Eq. ( |2.21| ) reads: 



Ai 



1 C-<pi 
2 v r 



1 

2Q 7 



A/Ajt 



A, 



(6.9) 



Purely self-dual Einstein vacuum spaces are obtained by demanding C- = (i.e. R fl f, = 
since the scalar curvature vanishes). This leads to the two known possibilities for Bianchi 



IX vacuum self-dual Einstein geometries met in Sec. 3.2 and satisfying either one of the 
following systems: 



Lagrange ( |3.13 l for A, = 0, 



Darboux-Halphen (3.14 i for A, = Cl l . 



6.2 Solving I & II with Painleve VI 



Systems I and II (Eqs. 1 6.7 1 and ( |6.8| ) describing general conformally self-dual Bianchi IX 
foliations with vanishing scalar curvature (W~ 



s = 0) were studied e.g. in |69 70J prior to 
their uplift to quaternionic spaces. Further developments in relation with modular proper- 
ties can be found in | [7T|[72| . 

As usual it is convenient to move to the complex plane, introduce o^(z) and 5^{z) and 



trade the dot for a prime as derivative with respect to z in (6.7 1 and ( |6.8| . Real solutions are 
recovered as previously: Cr(T) = ico { (iT) and A^(T) = iSi(iT). 

The system I is that of Darboux-Halphen for 5j(z) (see ( 3.14| ). Given a solution <5,(z) one 
can solve the system II for co'(z). Furthermore, the SL(2,C) solution-generating technique 



described in (4.2 1 can be generalized in the present case: given a solution Si(z) and a/(z), 



Si(z) 



(cz + d) 



2 Si 



az 



cz ■ 



+ 



cz 



(6.10) 
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and 



CO' z 



r£0' 



az + b 
cz + d 



(6.11) 



(cz + dy 

provide another solution if ( a h d ) G SL(2,C). 

Assuming S\ ^ 82 ^ 83 i.e. the triaxial situation (implying automatically to 1 7^ to 1 ^ co 3 ), 
we can readily obtain the general solution of the system I as in l|4.7}, 



Si(z) 



1 d_ 
2dz 



log5''( Z ), 



(6.12) 



with£ ! (z) a triplet of weight-two modular forms of T (2) C SL(2,Z). These can be expressed 
as in (J4T2I, where A is a solution of Schwartz's equation ( 4.11) . Define now a new set of 
functions Wi(z) as 

(6.13) 



IV; 



k cyclic permutation of 1, 2, 3), and insert the solutions (6.12 1 in system II ( |6.8) . The latter 
becomes 

dw\ W2W3 &W2 W3W1 dw?, W\W2 



dA 



A 



dA 



A-l' 



dA 



(6.14) 



A(A-l)- 

Notice the first integral w\ — w\ + a; 2 . Even though the value of this integral is arbitrary the 
uplift of the corresponding conformally self-dual geometry with zero scalar curvature to an 
Einstein manifold is possble only if the constant is V4 (see | |65}[66| ). 

The system of equations ( |6.14 1 can be solved in full generality with W[ expressed in terms 
of solutions y(A) of Painleve VI equation |73| (see also |74) for a more general overview): 



cl'T 



(y-A)y 2 (y-i) 



Here, 





A 


(y 


-A)y(y-i) 2 




A-l 


(y 


-A) 2 y(y-i) 




A(A-l) 




A(A-l)y' 



2(y-i) 
1 

~2y~ 
1 

2y 



2(y-A) 



2(y-A) 
1 

2^1) 



+ T- + 



2y(y-l)(y-A) 4y 4(y-l) 
and y is a solution of Painleve VI equation (/' = d //dA): 



4(3/ -A) 



(6.15) 
(6.16) 
(6.17) 

(6.18) 



, (y-A)y(y- 

^ 8A 2 (A-1) 2 



y- A 
l) 2 



(y 



A 



'\2 



1 1 

+ 



A A-l y-A 
- 1 3A(A - 1 



(y-i) 1 (y-A) 2 



(6.19) 
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6.3 Back to quasimodular forms 



We will for concreteness concentrate on the original solution of system I, the Halphen solu- 
tion corresponding to Ah = *f/ This is sufficient as any other can be generated by SL(2, C) 
transformations. Equations ( 6.14| (system II) read now 



w[ = irc&\w2W3, w' 2 = —indjio^wi, 



-ind\w2W2,. 



(6.20) 



In this form, the system can be solved in terms of Jacobi theta functions with characteristics 
|72[, as an alternative to the solution ( 6.15| -( 6Tl7| l. This makes it relevant in the present 
framework. 

2 — 1/4- required for the subsequent promotion to quater- 



The solution with w \ - 
nionic geometries - read: 



zvt; 



wi(z) 
w 2 {z) 

07 3 (z) 



27ti? 2 (0|z)j? 3 (0|z) 




(0|z) 






-J(0|z) 


27ri? 3 (0|z)i? 4 (0|z) 


m 


(0|z) 


_ e -i7t*/2 


9^CJ](0| Z ) 



27n? 2 (0|z)# 4 (0|z) m(Q\z) 



(6.21) 
(6.22) 
(6.23) 



Here a, b € C are moduli, mapped under the SL(2,C) transformations ( 6.10} and (6.11 



to 



other complex numbers. If a, b are integers and the transformation is in SL(2, Z), the solution 
is left invariant, up to permutation of the three components. 

It would be interesting to present the geometrical structure of the conformally self-dual 
zero-curvature spaces obtained with the solutions at hand, following the general procedure 



used in Sec. 5.2 This would definitely bring us far from the original goal. The interested 
reader can find useful information in the already quoted literature, both for these spaces and 
for their quaternionic uplift. Note in that respect that even though many families of solutions 
exist (here in the triaxial case, or more generally for biaxial three-sphere foliations), very 
few are singularity-free among which, the Fubini-Study or the Pedersen instanton (SU(2) x 
U(l) isometry), or the Hitchin-Tod solution (strict SU(2) symmetry). 

As a final remark, let us mention that ( 6.21| , (6.22 1, ( 6.23| also capture the self-dual Ricci 
flat solutions discussed in Sec. |3]and given in Eqs. ( 4.21| i.e. the Atiyah-Hitchin gravitational 
instanton. They correspond to the choice a = b = 1 mod 2, that must be implemented with 
care: consider a = 1 + 2e, b = 1 + 2zqb and take the limit e — > 0. One finds (a useful identity 
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for this computation is given in ( A.5 1: 



«H = ^o\^—--y [E 2 -4-n )), (6-24) 





( i 


TV 




\Z + Zq 


~6 




( i 


n 






~6 




( i 
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~6 



"4 = --kz\^-T[E2 + ei + Oi)}, (6-25) 



A modulus Zo is left in the solution; under SL(2, Z) it transforms as 

dzp + b 

z -> ■ — . (6.27) 

CZo + fl 

For finite zo the corresponding metric is Weyl-self-dual with zero scalar curvature The zq — > 
zoo limit corresponds to the Ricci-flat, Atiyah-Hitchin instanton (Riemann-self-dual). 

6.4 Beyond Bianchi IX foliations 

We would like to close our overview on conformally self-dual geometries with another fam- 
ily of quaternionic solutions, related to modular forms but not of the type M.± = R x M.^ 
with homogeneous M.$. Indeed, self-duality (Eq. ( 2.23) ) or conformal self-duality (Eqs. 



(2.24 1 and ( 2.25) ) can be demanded outside ot the framework of foliations. 



On can indeed assume an ansatz for the metric of the Gibbons-Hawking type fll"4) : 

ds 2 = O" 1 (dr + CDidx 1 ) 2 + <Mydx W. (6.28) 



Here O and CO, depend on x only, and thus 3 T is Killing. With this ansatz more general self- 
dual solutions are obtained with 11(1), 11(1) x U(l)orlT(l) x Bianchi isometry. Determining 
quaternionic spaces, i.e. conformally self-dual and Einstein, is however far more difficult. It 
is a real tour deforce to find the most general quaternionic solution with U(l) x U(l) isometry 



and this was achieved by Calderbank and Pedersen in [21 1, following the original method of 
Lebrun |68| . This will be our last example, where a new kind of modular forms emerge. 
In coordinates {p, rj,6,ip} with frame 

oc = VP~dp, p= gi+ggg 7 = dp/ S = dr j r (6 .29) 

VP 
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The metric reads: 



4p 2 ( F 2 + F 2 )-F 2 
= ^ P 4f2 / (7 2 + ^ 2 ) 

[(F - 2pF p ) a - 2pF^f + [(F + 2pF p ) /} - 2pF,, 



V - F 2 



F 2 



Here F p = d p F and F« = 9„F, where F(p,f]) is a solution of 



p 2 ( d 2 +d 2 \ F = 3 f _ 



(6.30) 



(6.31) 



The metric (6.301 has generically two Killing vectors, 9#, 9^, and F{p,rj) is a harmonic 
function on H2 with eigenvalue 3 / 4. Indeed, the metric on the hyperbolic plane is 



ds 



dp 2 + d?7 2 



H, 



(6.32) 



and Eq. ( |6.31 1 can be recast as 



A H2 F = f F. 



(6.33) 



Solving (6.331 leads inevitably to modular forms of r = y\ + zp, even though algebraic solu- 
tions are also available. 

Let us mention for example 



'P + 



(6.34) 



which leads to a metric on CP2 with U(l) x SU(2) isometry |18[, or 



P 2 -Pl 
2 VP 



(6.35) 



with U(l) x HeisenbergF] symmetry \75\. Solutions for F(p,rj) with strict U(l) x LT(1) isom- 
etry open Pandora's box for non-holomorphic Eisenstein series such as (see ( |A.9 |) 



F = £ 3 / 2 (t,t), 



(6.36) 



which has a further discrete residual symmetry SL(2, Z) C SL(2, R). These will be dis- 
cussed in Sec. [7] and we refer to the appendix for some precise definitions. Very little is 
known at present on the geometrical properties of the corresponding quaternionic spaces, 
or on the fields of application these spaces could find in physics. In string theory, they are 
known to describe the moduli space of hypermultiplets in compactifications on Calabi-Yau 



9 Heisenberg algebra is Bianchi II. 
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threefolds (76). The relevance of the Calderbank and Pedersen metrics in this context was 
recognized in (77). For further considerations on the role of modular and quasimodular foms 
as string instantonic contributions to the moduli spaces of these compactifications, we refer 
to 1 78 -83) and in particular to the recent review [84 1. 



7 Beyond the world of Eisenstein series 

To end up this review we would like to elaborate on some connections between quantum 
field theory and modular forms. This originates from the specific structure of the pertur- 
bative expansions in string and field theory, and calls for develping more general modular 
functions than the non-holomorphic Eisenstein series discussed earlier in these notes. 

7.1 The starting point: perturbation theory 

Perturbative expansions in quantum field theory are expressed as sums of multidimensional 
integrals, obtained by applying Feynman rules or unitarity constraints. These integrals are 



plagued by various divergences that need to be regulated. It was remarked in [85.86] that 
the coefficients of these divergences are given by multiple zeta values in four dimensions. 
Since this original work, it has become more and more important to further investigate the 
relationship between quantum field theory and the structure of multiple zeta values. This 
connection has fostered important mathematical results as in instance [87,88], which have 
been reviewed in the recent Seminaire Bourbaki by Pierre Deligne [89|. 

The next observation is that the above mentioned field-theory Feynman integrals arise 
as certain limits of string-theory integrals defined on higher-genus Riemann surfaces. They 
are actually obtained from the boundary of the moduli of higher-genus punctured Riemann 
surfaces. This bridge to string theory sets a handle to the world of modular functions. 

There are indeed two motivations for embedding the analysis into a string theory frame- 
work. The first is of physical nature: perturbative string theory is free of ultraviolet diver- 
gences, so it provides a well-defined prescription for regularizing field-theory divergences. 
In other words, string theory acts as a specific regularization from which we expect to learn 
more on the fundamental structure of the quantum field theories. The second motivation 
is directly related to the topic of this text: string theory is the ideal arena for exploring the 
number theoretic considerations of quantum field theory and their close connection with 
modular forms. 

In the present notes we will focus on the case of the tree level and genus one, following 
the string analysis in [22,90] and the mathematical analysis in |9l) . We will explain in partic- 
ular that (non-holomorphic) Eisenstein series are not enough for capturing all available in- 
formation carried by the integrals under consideration. The presentation will be schematic, 
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aiming at conveying a message rather that providing all technical details. For the latter, the 
interested reader is referred to the quoted literature. 

Let us consider the following integral defined on the moduli space of the genus-g Rie- 
mann surface with four marked points: 



/dp/ fld^expf £ 2oc'k i -kjP(z i ,z j )), (7.1) 

JM g JT, g i=1 \l<i</<4 / 



where Aig is the moduli space of the closed Riemann surface Eg of genus g. There are four 
punctures whose positions z ; - are integrated over. We have introduced k\ + ki + ko, + k^ = 
with k[ ■ kj = representing external massless momenta flowing into each puncture. We 
will also set the Mandelstam variables s = 7k\ ■ k^ = 1k-$ ■ k^, t = 2k\ ■ k^ = Iki ■ kj, and 
u = 2k\ ■ ko, = Iki ' ki, obeying s + 1 + u = for the massless states at hand. The physical 
scale is the inverse tension of the string a'. 

The propagator or Green's function P(z, w) is defined on this Riemann surface by 

= jf d 2 Zy/=gP(z f w), (7.2) 

d z d 2 P(z,w) = 2nS^(z)- 2 ! C 2 8z " (7.3) 

g 

d z d w P(z,w) = -27iS^{z) + 7zJ^cv I {z)(^mCl)Y } 1 cv ] (w) / (7.4) 



where the ds 2 = g zz dzdz is the metric on the Riemann, Q the period matrix, and coj with 
1 < I < g the first Abelian differentials. 

7.2 Genus zero: the Eisenstein series 

At genus 0, i.e. for the Riemann sphere, the propagator is simply given by 

pW(z,w) = log|z-w| 2 , (7.5) 



and the integral in (7.1 1 can be evaluated to give 



i r (i + x's)t (i + «'t)r (i + x'u) 



tx ,3 stu r (i - u's) r (l - cc't) r (i - x'u) 



n=\ 



The masses of string theory excitations are integer, quantized in units of !/«'■ It is therefore 
expected that the a! expansion of the string amplitude in ( |7.6| l is given by multiple sums over 
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the integers, but it is remarkable that this expansion involves only odd zeta values of depth 
one. For «' <C 1 a series expansion representation reads: 



A^\s,t,u)= £ c {M) a v 2 al, (7.7) 
?>-i,p>o 

where we have introduced cr 2 = (a's) 2 + (oc't ) 2 + (ol'u) 2 and 03 = 3a /3 sfu. Since <T\ = (a's) + 
(a'f) + {pc'u) = 0, we immediately see that all a n = (cc's) n + (a!t) n + (cc'u) n with n > 2 are 
given by 1 90 1 

£ = E il ^(f) P (f)'- <™> 

The coefficients cv p „) are polynomial in odd zeta values of weight 2p + 3q — 3. It is no- 
table that at a given order n = 2p + 3q — 3, the space of these coefficients has dimension 
d n = [{n + 2)/2j — [(n + 2)/3j, which coincides with the dimension of the space of the 
holomorphic Eisenstein series of weight n (see appendix). This hint calls for further inves- 
tigation, and we would like to mention the recent work connecting the a! expansion in p-7) 
and the motivic multiplet zeta values [92 J. 



One can expand the integrand of (7.7 1 and obtain each coefficient C( p „) as a linear combi- 



nation of the multiple integrals of the propagator p(°) (z,, Zy) (given in ( 7.5 ): 



'f!l2/"l3/"l4/"23/"24/«34 



/ n ^ n p {o) (^ Zj r>. (7.9) 

2 l<i<;<4 i<!<;'<4 



The integrand of this expression is the product of the propagators connecting the punctures 
with multiplicities n,y, 1 </<;'< 4, as depicted on Fig.|5j The contributions in ( |7.7[ | are the 
lowest-order to the full string-theory amplitude of the four-point (four punctures) processes 
described here. 



In Eqs. (|7.6[H 7.7 1, we encountered the zeta values 



W>1 " 

Extending the sum over the integers n to a lattice like p = m + rn G A^ 1 ' = Z + tZ, where^j 
t 6 [) = {z 6 C, 9m(z) > 0}, one gets the (non-holomorphic) Eisenstein series 

Ux,r)= E * (7.11) 

rZ IP I 



where |p| 2 = (m + nr){m + nr) is the natural Euclidean norm on the lattice A^. This 



10 The modular parameter T is expressed in alternative ways throughout these notes: T = 5Re(r) + i^sm(r) 
t x + ix 2 = rj + ip. 
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Figure 5: Graph of a vacuum Feynman diagram on the Riemann surface Eg. The punctures 
are connected by nu > links representing the number of two-dimensional propagators. 



expression can be made modular-invariant in a trivial way by multiplying by $5Wz(t) s , 

3?«(t) s 



E s (r / f)= £ 



|2s 



(7.12) 



This Eisenstein series is an eigenfunction of the hyperbolic Laplacian (6.31 1 with eigenvalue 

s(s-l): 

A H2 E s (T / f)=s(s-l)E s (T / f). (7.13) 



The case s = 3 /2 was discussed in Sec. 6.4 Eqs. ( 6.34| >— (6.36 l r from a different physical per- 
spective. 



7.3 Genus one: beyond 

At this stage of the exposition the reader may wonder how the above generalization of the 
zeta values ( 7.10} into modular forms ( 7.12} arises in string theory. We will sketch how this 
goes and show that new automorphic forms are actually needed, standing beyond the well 
known Eisenstein series. This requires going beyond the sphere (7.5l-( [7T7| |. 
At genus one, the Green's function is given by 



P«( Z ,0) = --log 



0i(z|t) 



^(0|t) 



2t 2 



(7.14) 
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The amplitude in ( |7.1| reads: 

A^(s,t,u)= [ ^ [ FT ^wWe-^^^W' 1 '^ (7.15) 

^ T 2 ^l<!</<4 T 2 

where J 7 = {r; |Ke(r)| < 2>^ m ( T ) > 0,3?e(T) 2 + Q'mfV) 2 > 1} is a fundamental domain 
for SL(2,Z), and T = {z;|Ke(z)| < |,0 < 3m(z) < 9m(t)}. 

No closed form for the integral in ( |7.15| l is known, in particular because of the presence of 
non-analytic contributions in the complex (s, £)-plane. For a rigorous definition of this inte- 
gral we refer to |93) . The expression for the genus-one propagator in ( 7.14| | has an alternative 
representation: 

pW(z,0) = — X] rT2 e_7r ~ +C(T,f), (7.16) 

Z7r peZ+rZ IPI 

where C(t, t) = log | ylrctj (t) | is a modular anomaly. Since the latter is z-independent, it 



drops out of the sum in (7.15 because of the momentum-conservation condition Yh=\ %i = 0- 



The integrand of (7.151 is therefore modular-invariant. From now on we will only consider 



the modular-invariant part of the propagator 

? (1) M) = ^ E ^e-^. (7.17) 
Following the developments on the sphere, we can analyze the expansion of the ampli- 



tude (7.151 for a.' <C 1. In this regime one gets integrals of the type §J B\ , but this time with 



the genus-one propagator 

Dn 12 ,n mW34 (T^)=/ ]1 Ft P {l KzuZj) n,i - (7-18) 

' l<«</<4 T2 l<i</<4 

The product runs over the entire set of links with multiplicities nu, 1 < i < j < 4, of the 
graph T depicted in Fig|5j By construction these integrals are modular functions for SL(2, Z). 
Performing the integration over the position of the punctures one gets an alternative form 
for the modular function D„ 12 . v „ 34 (t, f ) given by 

Dn 12 n»(T,T) = E Yl*( E Pi] EI (7-19) 



P,er i=i Xj—tvi / proper 



IP/ 



where the sum is over all the propagators p, of the graph T. If there are nyi propagators 
connecting the vertices 1 and 2, we have n\2 different elements of the lattice p, = m, + TWj G 
Z + tZ, 1 < i < n\i- At each vertex V\ of the graph we impose momentum conservation by 
demanding that the sum of the incoming momenta py flowing to this vertex (j — > vf) be zero. 
This is represented by the delta function constraint S(^j_y v . py) with 5{m + rn) = S(m)S(n). 
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The above sums D ni2( ..„ n34 (t, f ), introduced in | |22) , are generalizations of the Eisenstein 
series, that we will call Kronecker-Eisenstein following [91]. With each modular function 
D f!l2/ ... /f!34 (t, f ) we associate a weight given by the sum of the integer-valued indices n,;. Let 
us focus for concreteness on the particular case of n propagators between two punctures, 
and refer to p2) for the general case. We define 

D «(^):= E *(X>)nSf (7-20) 

P;£Z + TZ V ! ' =1 / ! = 1 

The special cases n = 2 and 3 are giverp]in [22 , appendix B] 



D2(T ' f) = w 1 ' (7 ' 21) 

However, in general these modular functions do not reduce to Eisenstein series, as it can 
easily be seen by evaluating the constant terms. For n > 2 it is always possible to decompose 
the modular form D„(t, f ) as |22j 

D n (t, f ) = P n (E s (t, t) ) + 5 n (t, t) (7.23) 



with P„(E s (t, t)) a polynomial in the Eisenstein series E s (t, t) (see Eqs. ( |7.12| and ( A.9| l) of 
the form 

P„(E 8 (T / f)) = p„(C(2n + l)) + ^^E J ,(T / f)+ E j^E r (r,f)E s (r,t), (7.24) 

where p n (^(2n + 1)) is polynomial of degree two in the odd zeta values of total weight n. 
The remainder 5 n (r, f) in ( |7.23[ ) is a modular form whose constant Fourier coefficient does 
not vanish but tends to zero for T2 — > 00. 

Although the definition of the modular functions D„ 12/ ... /(!34 (t, f ) given in §7.20) looks sim- 



ilar to the double-Eisenstein series introduced in [94|, one finds that, as opposed to the latter, 
their constant term involves depth-one zeta values [22[ only. Hence, they provide a natu- 
ral modular-invariant generalization of the polynomials in the odd zeta values met in ( |7.9| . 
One way to obtain multiple zeta values is to insert in ( |7.20| | the generalized propagator used 
by Goncharov in |9l| . Whether the generalization introduced by Goncharov does appear 
in string theory is an open question. From the original physical perspective, this question is 
relevant because it translates into the (im)possibility of appearance of multiple zeta values as 
counter-terms to ultraviolet divergences in quantum field theory. This might have important 



11 The case n = 3 has been worked out by Don Zagier. 
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consequences in supergravity. 

As a final comment, let us mention that although our discussion was confined to the case 
of modular functions for SL(2,Z), most of the above can be generalized to the framework 
of automorphic functions for higher-rank group |23|. 



8 Concluding remarks 

In the present lecture notes we have given a partial - in all possible senses - review of the 
emergence of (quasi)modular forms in the context of gravitational instantons and string the- 
ory. These forms often appear as the consequence of remarkable, explicit or hidden symme- 
tries, and turn out to be valuable tools for unravelling a great deal of properties in a variety 
of physical set ups. The latter include monopole scattering, Ricci flows, non-perturbative (in- 
stantonic) corrections to string moduli spaces (via their Fourier coefficients), or perturbative 
expansions in quantum field theory (via string amplitudes). 

We have described how the classical holomorphic Eisenstein series, whose theory is 



nicely presented in [95 .96], occurs in the context of gravitational instantons or in studying 



non-perturbative effects. We have also encountered the non-holomorphic Eisenstein series, 



the analytic properties of which are described in [97,98]. This whole analysis has led us to 
argue that one needs novel types of modular functions, standing beyond the usual Eisen- 
stein series. Although the analytic properties of these series are still poorly understood, they 
seem to be a corner stone for understanding the challenging nature of interactions in string 
theory and its consequences in quantum field theory. 
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A Theta functions and Eisenstein series 



We collect here some conventions for the modular forms and theta functions used in the 



main text. General results and properties of these objects can be found in [95 96 [. 
Introducing q = exp 2inz, we first define 



n=l 
12 d 



fj{z) 

E 2 (z) = "-^ log 77, 
ik dz 



as the Dedekind function and the weight-two quasimodular form, whereas 



(A.1) 
(A.2) 

(A.3) 



are the Jacobi theta functions. More generally, one introduces 



(v\z) = ex P [}Tiz{m + «/2) 2 + 2in(v + b /i){m + a /i)) 



(A.4) 



with ti[\] = = $2,$lo] = % #[?] = #4- Let us also mention the following relation 



a. + 2w 



(v\z) 



(v + wz\z). 



(A.5) 



The first holomorphic Eisenstein series are 



E 2 (z) = l-24E^ 1 ^ 
E 4 (z) = l+240E£ =1 ^ 



3,,m 



-q rn 

2*(z) = l-504nUg£. 



(A.6) 



Notice that E4 (z) and Eg (z) are modular forms of weight 4 and 6, whereas E2 (z) is the already 
quoted weight-two quasimodular form. The modular-invariant of weight two is the non- 
holomorphic combination £2(z) — 3 / 'nQm(z). It is a classical result that the space of modular 
forms of weight k is spanned by E^Eg with 2a + 3b = k. The dimension of this space is 
d fc =L(* + 2)/2]-L(* + 2)/3j. 

In the main text we also consider non-holomorphic Eisenstein series E s (z, z) with z = 
x + iy, y > and x S R. These are defined as modular-invariant eigenfunctions of the 
hyperbolic Laplacian (see Eqs. ( 6.31) , ( 7.12) and ( 7.13) ) 



y 2 {d 2 x + 3J) E s (z,z) = s(s - l)E s (z,z), 



(A.7) 
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with polynomial growth at the cusps (y — > oo): 



e s (z,z)= £ i — ^n*,* < A - 8 > 
(m,nmm \ mz + n \ 2s 

for s G C with large enough real part for convergence. One can extend the definition by 
analytic continuation [97] for all s ^ 1 using the functional equation T(s)n~ s E s (z,z) = 
r(i/2 - s)n 1/2 - s Ei_ s (z,z). These series have the following Fourier expansion: 

2£(2s)E s (z,£) =2^2s)y s + 2£(2s-l)i/ 1 - s +4i / 1 / 2 £ P "| S "^ 2 l) X s _ 1/2 (27i|n|y) e 2 ^ (A.9) 

where £(s) = £,(s)T( s /2)n~ s/2 is the completed zeta function, K s _i/ 2 is the K-Bessel function 
and c a (n) = YLd\n d a (see e.g. |99| for details). 

Finally let us mention how the non-holomorphic series are connected to the holomorphic 
ones. For that, one considers the following generalization of the non-holomorphic Eisenstein 
functions: 

..s+w+w 

E ( f' w \z,z)= Y -. L (A .10) 

(mz + n) s+w {mz + n) s + w v ' 



s 



(m,n) ^(0,0) 

c d , 

: {w,w), \ _ ( j-sw i j,w r (w,iv) , 



These series transform under a modular transformation 7 = f a \ ) G SL (2, Z) as 



'(7-z,7-f) = (cz + d) w (cz + d) a; Er' J (z,z). (A.ll) 

Chosings = n G Nando) = —n, we recover the holomorphic Eisenstein series E,^, ' (z,z) = 
E„(z). 
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